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Abstract 

At present, approximately 1500 asteroids are known to evolve inside or sticked 
to the exterior 1:2 resonance with Mars at a ^ 2.418 AU, being (142) Polana 
the largest member of this group. The effect of the forced secular modes 
superposed to the resonance gives rise to a complex dynamical evolution. 
Chaotic diffusion, collisions, close encounters with massive asteroids and 
mainly orbital migration due to the Yarkovsky effect generate continuous 
captures to and losses from the resonance, with a fraction of asteroids re- 
maining captured over long time scales and generating a concentration in 
the semimajor axis distribution that exceeds by 20% the population of back- 
ground asteroids. The Yarkovsky effect induces different dynamics according 
to the asteroid size, producing an excess of small asteroids inside the reso- 
nance. The evolution in the resonance generates a signature on the orbits, 
mainly in eccentricity, that depends on the time the asteroid remains cap- 
tured inside the resonance and on the magnitude of the Yarkovsky effect. 
The greater the asteroids, the larger the time they remain captured in the 
resonance, allowing greater diffusion in eccentricity and inclination. The res- 
onance generates a discontinuity and mixing in the space of proper elements 
producing misidentification of dynamical family members, mainly for Vesta 
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and Nysa-Polana families. The half-life of resonant asteroids large enough 
for not being affected by the Yarkovsky effect is about 1 Gyr. From the point 
of view of taxonomic classes, the resonant population does not differ from 
the background population and the excess of small asteroids is confirmed. 

Keywords: Asteroids, dynamics. Asteroid Vesta, Resonances, orbital. 
Rotational dynamics 



1. Introduction 

The asteroid Main Belt is sculpted by mean motion resonances (MMR), 
mainly with Jupiter, and secular resonances with the giant planets. However, 
there are also relatively strong exterior MMR with the terrestrial planets. 



especially in the inner region of the Belt (jGallardd . l2006l ) . For any terrestrial 
planet, its exterior 1:2 MMR is the strongest one, and in the case of Mars, 
this resonance (1:2M) is located in the inner main belt at a ~ 2.4184 AU. 
This MMR prod uces an evident concentration of asteroids in semimajor axis 



(jGallardd . 120071 ). How was this concentration created? How long does an 
asteroid stay in resonant motion? Is there any mechanism that replenishes 
the population? The study of the dynamics of the 1:2M resonance and the 
origin and evolutionary paths of their members is the aim of this paper. 

If we try to approximate the perturbing function for this resonance by 
an analytical approach using a series of terms of different orders in the small 
parameters e and i, the lowest order terms depend on the critical angles a = 
2X — X]\f — tu and ai = 2X — Xm — ^m, where subscript m refers to Mars. These 
terms are the ones that are usually librating in the case of asteroids located 
in the resonance. The first one is factorized by the asteroid's eccentricity, 
e, and the second one is factorized by Mars' eccentricity, Cm- Higher order 
terms, and consequently dynamically less important ones, depend on other 
combinations of A, Q, Xm, Wm and VLm-, which define other critical angles. 
In very particular cases, these critical angles can librate too. 

In this resonance, as in any exterior resonances of the type 1:N, the libra- 
tions of the critical angle a occur around central values, dc, that depend on 
the eccentricity, giving rise to the so-called "asymmetric librations" because 
they do not occur around 0° or 180° as is usual in other MMRs. The angles 
a and cii, and especially the first one, are the most important indicators of 
the resonant state but not the only ones. In particular, the libration theory 
predicts an oscillation of the osculating semimajor axis around a mean, a^, 
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or "proper" value, Cp, which is characteristic of the resonance. The difference 
between and is that the first one is the mean over some time span and 
the latter is obtained from a frequency analysis. Moreover, it is known that 
the resonance width, i.e. the region in semimajor axis where the resonant 
motion dominates, is proportional to the orbital eccentricity defining a region 
in the (a, e) space where the resonance drives the orbital evolution. 

Approximately 1500 asteroids are performing librations, horseshoes and 
transitions between both asymmetric libration centers, which are the typical 
trajectories in this resonance. About 400 of these asteroids librate around 
the asymmetric libration centers with amplitudes Act = a^ax — o'mm < 180°. 
The states of the resonant populati on now a n d afte r 1 Myr of orbital evolu- 
tion are shown in Figs. 1 and 2 of iGallardd (120091 ). Asteroids change their 
libration center and change between librations and horseshoe trajectories, 
but the number of asteroids with asymmetric librations and horseshoes is 
similar at the beginning and at the end of this time-span. Therefore, we 
may conclude that the population is in equilibrium over such timescale. Ta- 
ble [T] lists the resonant asteroids with absolute magnitude H < 14. The 
evolution of the critical angle and eccentricity of (142) Polana is shown in 
Fig. [U Proper elements shown in the table were derived from our numer- 
ical study which is explained in Section [21 and the complete list can be 



found at http://www.fisica.edu.uy/~gallardo/martel2/all.html In 



some cases, normally at low e, the libration amplitude of ai is smaller than 
the libration amplitude of a. 

At variance to the resonances with Jupiter that dominate the time evo- 
lution of the asteroid eccentricity, the libration theory predicts that the ec- 
centricity is almost unaffected by the 1:2M resonance, but its evolution is 
dominated by the secular forced modes on the asteroid orbit generated by 
the planetary system. In fact, the variations in eccentricity shown in Fig. [T] 
are due to a secular evolution, not to the libration itself. These important 
oscillations in the eccentricity affect the locations of the libration center and 
the ampl i tude o f the asymmetric librations of the critical angle (Fig. 2 from 



Gallardd . 120071 ). Nevertheless, these are long term periodic oscillations and 



they do not destroy the resonant motion. This can be checked by follow- 
ing the time evolution of Cp or am, which remain locked in the resonance's 
domain. However, the irregular evolution of a no longer guarantees the res- 
onant protection mechanism against close encounters with Mars, allowing 
such close encounters for high eccentricity orbits. 

The region around the location of the resonance is occupied by the dynam- 
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Asteriod 


H 






hi] 




Aa[°] 


(142) Polana 


10.2 





1573 


3.198 


105 


211 


(1998) Titius 


12.2 





0848 


7.764 


70 


124 


(doooj ritzgerald 


12.6 





1151 


9.557 


289 


95 


(9652) 1996 AF2 


13.0 





1967 


6.860 


96 


125 


(2798) Vergilius 


13.1 





0438 


5.944 


214 


200 


(11576) 1994 CL 


13.3 





1394 


11.111 


80 


131 


(11055) Honduras 


13.5 





1937 


11.774 


70 


72 


(11751) 1999 NK37 


13.7 





1347 


2.656 


94 


142 


(8748) 1998 FV113 


13.7 





0917 


5.937 


226 


170 


(42786) 1998 WU4 


13.8 





1196 


22.383 


141 


197 


(2994) Flynn 


13.9 





1974 


2.649 


275 


102 



Table 1: The largest resonant asteroids with critical angles (a) showing asymmetric libra- 
tions (libration center (Jc 0°, 180°) or transitions between libration centers. H is the ab- 
solute magnitude, Cp and ip are proper eccentricity and inclination, and Act ~ (Tmax — Cmin • 



ical families of Vesta, Massalia and the Nysa-Pqlana complex (jCellino et al. 



20011 : iMothe Diniz et al.l . l2005l : IZappala et al.l . [19951). The 1:2M resonance 



gener ates some dynamical signatures to the Massalia family (jVokrouhlicky et al. 



20061 ). and should generate some dynamical signatures to the other families, 



a topic that we will explore in this paper. 

The resonance is very narrow in semimajor axis units, but even so it 
generates a feature in the distribution of main belt asteroids, which means 
that it is strong enough to affect their dynamical evolution. We are interested 
in determining the resonance ability to retain asteroids in spite of the chaotic 
diffusion, the Yarkovsky effect, and the mutual collisions and close encounters 
with massive asteroids that may occur in that region. 



2. Resonant population: dynamical and physical properties 

We have explored the dynamics of the res onance by a series of n umeri- 
cal integrations using the integrator EVORB ( iFernandez et al.l . 120021 ). Our 
model includes the planets from Venus to Neptune, ignoring the planet Mer- 
cury, plus real and fictitious populations of asteroids. The orbital elements 
were taken from the ASTORB database (ftp : // ftp . lowell . edu/pub/ elgb/ astorb . html ) 
and the output of the integratio ns was analyzed using fourier techniques 
(IGallardo and Ferraz-Melld . 119971 ). The numerical integrator is basically a 
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leapfrog scheeme with time step of 0.02 yr in the absence of close encounters, 
which switches to a Bulirsch-Stoer routine in the case of close encounters of 
the asteroids with the planets. We also analyzed the physical properties of 
the asteroids, like absolute magnitude, if, surface colors from the Sloan Dig- 
ital Sky Survey Moving Objects Catalogue (SDSS-MOC) photometry, and 
membership to the known dynamical families in this region. 

2.1. Resonance dynamics 

In order to find the properties of the resonant motion, we looked for 
known asteroids with osculating semimajor axes near the resonant value and 
follow their orbital evolution by means of numerical integrations. We started 
our analysis by integrating 4022 orbits from the ASTORB file with 2.416 < 
a < 2.422 AU. The resonance domain in semimajor axis can be visualized 
in several ways. For example, in Fig. [21 we show a plot of Cp versus the 
principal frequency associated with the time evolution of a. Both values 
were deduced automatically from the spectral analysis of the time evolution 
of a and a over 30 x 10^ years. We define as the time-independent term 
in the fourier series expansion of the osculating a{t). In Fig. |21 the central 
feature at Op ~ 2.4184 AU corresponds to asteroids performing librations 
with periods going from 2,500 to 10,000 yr. The two inclined lines at both 
sides of the central feature are due to asteroids with their critical angles 
performing circulations. Moving in Op from the left (where da/dt > 0) to the 
right, the circulation frequency diminishes up to the point where da/dt ~ 0, 
and the librations start generating the discontinuity in Cp which accumulates 
at Op ^ 2.4184. Continuing to the right we find, again, the circulation of the 
critical angle but in this case with da/dt < and with increasing frequency. 
Between the central feature and the inclined lines, there are two regions where 
transitions between circulation and libration are observed. For example, an 
asteroid with Op = 2.416 shows only circulation in the critical angle, but this 
circulation frequency is evident in the time evolution of a, which means that 
there are traces of the resonance in its dynamical evolution. No frequencies 
associated with a were found for Op < 2.414 and Op > 2.423, so no dynamical 
traces of the resonance can be expected there. 

Performing a spectral analysis of the variables [k = e cos to, h = e sin w) 
and [p = icosQ, q = isinQ) over 1 Myr of numerical integration, it was 
possible to determine the forced and proper modes in the evolution of e and 
i. Forced modes are the oscillation modes corresponding to the fundamental 
frequencies of the planetary system, and proper modes are the oscillation 
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modes wi th large amplitudes and fre quencies different from the fundamen- 
tal ones ( Murray and Dermott . 19991 ) . We call and ip the amplitudes of 



the forced modes. Contrary to Fig. [21 where we cannot distinguish be- 
tween different eccentricities, Fig. [3] shows the structure of the resonance 
with increasing width at higher Cp, and the huge concentration of asteroids 
at 2.4180 < flp < 2.4188, with two deep gaps at both sides of the cen- 
tral peak. An analogous plot using synthetic proper elements from AstDyS 
(Ihttp : //h amilton.dm.unipi . it/astdys/) looks very similar to our Fig. [3l 
Considering the peak and the gaps, there is an excess of asteroids in the 
region 2.416 < Op < 2.421. According to the AstDyS database, that only in- 
cludes asteroids with well determined orbits so it should be less biased than 
our sample, we estimate that in a range of Aa^ ~ 0.005 AU there are approx- 
imately 1900 asteroids in the domain of the resonance, against approximately 
1600 asteroids that are far from the resonance (background). Therefore, in 
that database, there is an excess of approximately 300 resonant asteroids, 
which means that the resonance generates a population excess of 20% over 
the background. 

Figures |5] and |2] can be used to define the limits where the resonance is 
dominant over the other periodic terms that affect the time evolution of the 
asteroid semimajor axis. However, the best indicator of these limits appears 
in Fig. m that was obtained by plotting the 1 Myr averages of the orbital 
elements of 100 resonant asteroids from the ASTORB file evolving over 1 Gyr. 
The concentration at the borders of the resonance is due to the phenomenon 



of "resonance stickiness" ( iMalishkin and Tremaind . Il999[ ). and the behavior 



at Cm > 0.2 is due to the overlap of reso nances that cause a strong chaotic 



diffusion (IMorbidelli and Nesvornyl . 1 19991 ) . 



The fundamental frequencies and Qq are clearly present in the secular 
evolution of the eccentricity, and /e appears in the inclination. The proper 
elements Cp and ip are concentrated around some values, indicating their 
relation to particular coUisional families (Fig. |5]). 

2.2. Physical properties 

According to the orbital evolution over 1 Myr, we defined a subsample of 
asteroids clearly captured in the resonance (2.4180 < ttp < 2.4188) and a sub- 
sample of asteroids clearly outside of both sides of the resonance {ap < 2.4154 
and ttp > 2.4214). In order to determine the fraction of these bodies that 
belong to asteroid dynamical families, we applie d the Hierarchical Cluster ing 



Method (following the procedure described in iMothe Diniz et al.l . |2005[ ) to 
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Family 


a,[AU] 


e 


p 


h\ 


1 


% of R 


% of NR 


Vesta 


2.25 - 


2.48 


0.08 - 


- 0.12 


5.7- 


7.5 


2.9 


6.2 


Nysa-Polana 


2.32 - 


2.47 


0.14 - 


- 0.20 


2.0 - 


3.2 


5.3 


9 


Massalia 


2.36 - 


2.46 


0.15 - 


- 0.17 


1.3 - 


1.7 


3.3 


3.7 



Table 2: Percentage of dynamical family members in the subsample of resonant asteroids 
(R) and in the subsample of non resonant asteroids (NR) . There is a deficiency of identified 
family members inside the resonance. 



the AstDyS catalog of analytical proper elements, and identified the major 
asteroid families in the neighborhood of the 1:2M resonance. We found that 
only 11.5% of the clearly resonant asteroids belong to the dynamical families 
of either Vesta, Nysa-Polana or Massalia. For the non resonant subsample, 
only 18.9% are members of any of the above families (see details in Table |2]). 
The families of Vesta and Nysa-Polana clearly have less members inside the 
resonance than outside. Re s ults f or Nysa-Polana family are consistent with 



Fig. 26 of iMorbidelli et al.l (jl995[ ). where a decrease in the density number 
of the family members is observed at a ~ 2.42. As we will show later, the 
lack of family members in the resonance is relate d to the diffusion in (e, i) 



that hinders the identification of family members (jVokrouhlicky et al.l . 12006 



Carruba and Michtchenkd . 120091 ). 



For both resonant and non resonant groups, we analyzed the distribution 
of asteroids' surface colors using the photometric data from Sloan Digital 
Sky Survey (SDSS). We did not find any remarkable differences between the 
groups, which led us to conclude that the resonance is not dominated by any 
particular asteroid taxonomic class. The analysis of the absolute magnitude 
(H) distribution confi rmed a tendency to higher H (smaller asteroids) for the 
resonant population ( IGallardol 120071 ). However, the resonant population is 
more eccentric, on average, than the non resonant one (Fig. [6]), and this might 
favor the discovery of smaller asteroids inside the resonance than outside. 
To correct this bias, we analyzed the H distribution of two samples of 855 
asteroids each, one from the resonant population and the other from the non- 
resonant population, having the same eccentricity distribution in the range 
< e < 0.25 (Fig. [7D- Using a Kolmogorov-Smirnov test for cumulative 
distributions (jPress et al.l . 119971 ). we found that both samples have different 
H distributions with a level of confidence of 98.7%, and that small asteroids 
are more abundant in the resonant sample. We will discuss this point later 
in Section [3l 
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2.3. Long term evolution of the resonant population 

We integrated numerically 100 real asteroids over 1 Gyr, and calculated 
the mean eccentricity and the mean semi-major axis over intervals of 1 Myr. 
Then, using the limits of the 1:2M resonance defined by Fig. HI we determined 
the number of asteroids captured in the resonance at each instant and the 
result is shown in Fig. [SI By fitting this data to an exponentially decaying 
function, we can deduce a half-life of ~ 2.5 Gyr. This value will be revised 
later in Section [31 

While evolving inside the resonance, the asteroids exhibit a diffusion 
mainly in 6^, but also in i^, which in general is roughly proportional to 
the lifetime in the resonance. A typical variation in is about 0.05 after 
some hundred million years of evolution inside the resonance, but some as- 
teroids may exhibit larger changes in e^, and variations up to 0.15 can be 
achieved for the asteroids that remain captured over 1 Gyr. We have also 
found that the larger proper inclination, the larger the diffusion in Cm and 
im- For example, for ip = 1.5°, typical of the Massalia family, the diffusion in 
eccentricity is 0.05 and in inclination is 1.2°, whereas for ip = 6.6°, typical of 
the Vesta family, the diffusion in eccentricity is 0.08 and in inclination is 1.6°. 
The success in the identification of resonant members of the Massalia family, 
in comparison to the members of the Nysa-Polana family, could be explained 
partially because of their low proper inclination and, in consequence, their 
lower diffusion in eccentricity and inclination. Figure M shows the evolution 
in {em,i'm) for the complete sample. When crosses the value ~ 0.22, the 
diffusion in am and im increases and the asteroids evolve in a chaotic region 
generated by the overlap of several resonances, eventually leaving the 1:2M 
resonance. 

In our simulation, the first asteroid had a close encounter with Mars after 
84 Myr, and after 200 Myr there were 5 asteroids encountering Mars, so the 
resonance is a clear source of Mars crossing asteroids. The first encounter of 
an asteroid with the Earth happened at approximately 100 Myr. Six asteroids 
in the simulation, most of them with initial i > 11°, ended by colliding with 
the Sun shortly after being captured into the 3:1 resonance with Jupiter (the 
most common situation), or in the 4:1 resonance with Jupiter, or evolving 
near the z/g secular resonance. 

Other MMRs very close to the 1:2M resonance are listed in Table [3] 
These resonances have much lower strengths, however, for higher values of e 
they grow in relevance and can compete with the 1:2M resonance, allowing 
a diffusion in a. We have also found captures in three body mean motion 
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2.4116 
2.4137 
2.4184 
2.4237 
2.4351 



Resonance a[AU] strength (e = 0.15) strength (e = 0.25) 
19Aj-6A 
-4:Xe + 15A 
— Ajvf + 2A 
22Aj - 7A 
-5\e + 19A 



4 X 10" 
6 X 10-6 
1 

2 X 10-^ 
1 X 10 
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1 X 10"^ 

5 X 10-3 

1 

1 X 10-^ 

6 X 10-^ 



Table 3: The strongest two body mean motion r esonances near a 2.418. The resonance 
strength was calculated following iGallardol (|2006l) and is given in units relative to the 1:2M 
resonance strength for two different values of the orbital eccentricity and for fixed z = 5° 
and ui = 60°. Note the increasing strength of the resonances relative to 1:2M for the more 
eccentric orbits. 



resonances ( iNesvorny and Morbidellil . Il998[ ). as for example the Sum + tin — 
6n = resonance at a ~ 2.4125, that involves the mean motions of Mars 
and Neptune. All these MMRs contribute to the chaotic diffusion observed 
at high eccentricities. 



3. Orbital migration by the Yarkovsky effect 

All the results presented up to this point were obtained in the framework 
of pure gravitational intera ction between the asteroids and the planets. It is 
known (IBottke et al.l . 120001 ) that the Yarkovsky effect generates a small extra 
acceleration due to the natural delay of thermal emission by the asteroid 
surface, and this effect could significantly modify the capture mechanisms, 
especially for asteroids with diameter < 10 km. 

We found that 80% of the 4022 asteroids located in the 1:2M resonance 
region have absolute magnitudes 14.9 < H < 17.8, which corresponds to 
diameters between ~ 0.6 — 6.0 km depending on the geometrical albedo 
considered. So, we are dealing with a population of asteroids for which the 
Yarkovsky effect cannot be neglected. This effect can be reproduced with a 
simplified rnodel b y adding a term to the asteroid's acceleration of the form 
fiRoig et al.l . l2007t ): 

where G is the gravitational constant, Mq the solar mass, a the osculating 
semimajor axis of the asteroid, r* and v are its heliocentric position and 
velocity, respectively, and a the assumed Yarkovsky drift in a that depends 
on several factors like diameter, rotational axis obliquity, rotational period. 
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and surface thermal conductivity, K. According to Bottke et al. (j2000u2006[ ). 
a runs from ~ 10"^ to ~ 10"^ AU/Myr for asteroids with radius from some 
meters to a few kilometers and values of K between 0.01 and 0.001. The 
most common variant is the Yarkovsky "diurnal" effect, that makes asteroids 
with prograde rotation to have d > 0, and those with retrograde rotation to 
have d < 0. In very particular cases, like asteroids with iron-rich surfaces 
and high values of the thermal conductivity, another variant known as the 
Yarkovsky "seasonal" effect (that we will not consider here) predominates 
and produces a secular negative drift in a (IBottke et al.l . l2006l ). 

We have performed a series of four simulations with four different assumed 
values of d: 10~^, 10""^, 10"^ and 10"^ AU/Myr (hereafter referred as simu- 
lations #1, #2, ^3 and ^4, respectively, see Table H]). In each simulation, 
we consider two fictitious population of 200 test particles initially located at 
both sides of the resonance and with values of (e, i) taken from the real pop- 
ulation. In order to force the particles to fall into the resonance, we assumed 
d > (i.e, prograde rotation) for the particles starting from "the left side" of 
the resonance, and d < (i.e., retrograde rotation) for the particles starting 
from "the right side". This allows us to detect any difference between left 
and right side populations regarding the capture into the resonance. In these 
simulations, we do not consider the abrupt and stochastic changes in d due to 
eventual reorientation of spin axes related to non catastrophic collisions. We 
also ignore the obliquity and rot ation rate evolution due to the YORP effect 
(jVokrouhlicky and Capekl . 120021 ). It must be noted that the simplified model 
used to simulate the Yarkovsky effect generates a spurious variation in the 
particles eccentricity, but this is negligible in magnitude. Indeed, if we cal- 
culate the total variation in e induced by the Eq, ([1]) during our simulations 
only by means of the Gauss planetary equations, we obtain ~ 10~^ — 10~^ 
which is small enough to not alter our conclusions. 

The simulations spanned a time scale long enough for the total drift in a 
to be larger than the resonance width (i.e., the asteroids would completely 
cross the resonance during the simulations provided that no captures occur). 
We also estimated the time scale of the collisional/gravitational lifetime, tejec, 
as the time necessary to eject an asteroid from the resonance due to colli- 
sions with other asteroids, or due to gravitational perturbations by close 
encounter with massive asteroids. To compute tcjcc, we simulated collisions 
using a Monte Carlo method. We a ssumed a typical power law population 



of projectiles (IFarinella et al.l . Il998[ ) with randomly oriented encounter ve- 



locities, taken from a Maxwellian distribution with mean 5.3 km/s. We 
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Simulation 


d [AU/Myr] 


R [m] 


AT [Myr] 


Aa [AU] 


tcjec [Myr] 


#1 


0.01 


1 


2 


0.02 


5 C 


#2 


0.001 


30 


20 


0.02 


10 C 


#3 


0.0001 


600 


200 


0.02 


90 C 


#4 


0.00001 


5000 


1000 


0.01 


1000 E 



Table 4: Parameters of four simulations using different Yarkovsky drifts a. Each simulation 
was composed of two fictitious populations of 200 bodies. Each population was initially 
located at e ach side of the reson ance. The typical radius of the bodies, R, was estimated 



according to iBottke et al.l (|2006f ) and assuming a thermal conductivity between 0.01 and 



0.001. Table gives the total time span, AT, of the simulations and the equivalent total 
drift, Aa. The ejection timescale, ^ejoo has a large uncertainty and it was estimated as 
the mean time necessary for the particles to be ejected from the resonance due to either 
gravitational perturbations by close encounters with massive asteroids (labeled as "E") , or 
collisions with other asteroids, including disruption events (labeled as "C"). 



also assumed an impulse 6v = (-Rproj/-Rtargct)^'ye5 and co mputed 6a via Gaus s 



planetary equations, considering disruption events as inlBottke et al.l (12005 ) 



Gravitational perturbations by massive asteroids (iNesvorny et al.l . |2002[ ) 
were also estimated via Gauss planetary equations, using a Monte Carlo 
method to calculate the impulse 6v = 2mGve{G'^{M + raf + p^f^)"^/^ due 
to hyperbolic encounters with the 10 most massive asteroids, where p is the 
randomly chosen impact parameter, M the mass of the target assuming a 
density of 2.7 g/cc, and m the mass of any of the randomly chosen massive 
asteroids. The mean time interval between encounters with the 10 most 
massive asteroids and with impact parameter p < 0.01 AU was taken to be 
15,600 yr. This time scale was deduced from the intrinsic collision probability 
of 2.86 X 10~^^ km~^yr~^, typical for the inner region of the Main Belt. 

The gravitational perturbations by massive asteroids are approximately 
the same for all the range of asteroid sizes. On the other hand, the effect 
of collisions is only relevant for smaller asteroids. In our simulations, we 
found that tejec is determined by collisions for asteroids with R < 2.5 km 
and by close encounters with the massive asteroids otherwise. These close 
encounters drop off the half-life obtained from Fig. |8]to ~ 1 Gyr. 

3.1. Capture in resonance and chaotic diffusion 

The analysis of the time evolution of the semimajor axis shows that the 
main mechanism that keeps the asteroid fixed in a is the resonance stick- 
iness. Figure [10] shows a typical example of a particle from the "left side" 
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group in simulation ^4. The Yarkovsky effect generates a positive drift in 
semimajor axis up to t ^ 500 Myr, when the particle gets sticked to the res- 
onance alternating between the "left" and "right" sides. At t ~ 750 Myr, the 
particle gets captured into the resonance, and after that it evolves switching 
between sticking and getting inside the resonance. It is worth noting that 
after t ~ 500 Myr the drift in a stopped, and a new dynamical mechanism 
drives the eccentricity. This mechanism partially erases the "memory" of the 
particle, mainly its proper eccentricity but also its proper inclination. Cap- 
tures from the "right side" exhibit a similar behavior. In some cases, the 
particle is neither captured into the resonance nor in a sticky orbit, and in 
these cases the eccentricity is not altered. This is illustrated from the results 
of simulation #4 shown in Figs. [TT] and [T^l Figure [TT] shows the popula- 
tion evolving from the "right side" in the [am, Gm) plane, and Fig. [12] shows 
the population evolving from the "left side" in the {emiim) plane; the diffu- 
sion in this later case is significant for asteroids captured into resonance and 
negligible for asteroids that quickly escape from the resonance. 

Calculating the of each particle in a 1 Myr running window, we com- 
puted the number of asteroids inside the resonance or sticked to the resonance 
as a function of time. We also computed the time t^^s that each asteroid re- 
mains captured into/sticked to the resonance. We considered the asteroid was 
captured/sticked when falls inside the resonance region in the (am; e^) 
space defined in Fig. |H Analyzing all four simulations, we found that the 
total effect on the eccentricity when evolving inside the resonance is somehow 
proportional to the magnitude of the Yarkovsky effect and to the time the 
particle remained captured into the resonance (Fig. [T3l) . The diffusion in 
mean eccentricity, Acm, as a function of the lifetime in the resonance tres, in 
Myr, can be approximately described by 

Ae„, ~ 6 d (2) 

where a is in AU/Myr and /3 is 0.54, 0.75 and 1.07 for simulations ^2, 7^3 
and #4, respectively. The value of Ae,m deduced from this expression is only 
an indication of the actual value, considering the large spread of the data 
observed in Fig. [131 From this figure, we can conclude that, in general, 
the largest diffusion in eccentricity corresponds to simulation ^4, i.e., to the 
largest asteroids. 



Simplified analytical models (see for example [Murray and Dermottl . 119991 ) 



predict that captures from "approaching" orbits (i.e., asteroids migrating 



12 



towards Mars, or equivalently, from the "right side" population) are more 
successful than captures from "receding" orbits (i.e., from the "left side" pop- 
ulation), but in our simulations we did not find any significant difference 
between the two populations regarding the probability of capture. As an 
example. Fig. [H] shows the fraction of asteroids captured into resonance as a 
function of the time in simulation ^3 (i.e., for a Yarkovsky effect correspond- 
ing to an asteroid size of some hundred meters), distinguishing between the 
left and right populations. The same behavior was obtained in all the four 
simulations, but in different time scales. 

Using the values of tres computed for all particles in all the four simula- 
tions, we determined the cumulative distribution of the scaled lifetime in the 
resonance, L = tj-^s/T, where T is the time that would be necessary to cross 
the resonance region due to the Yarkovsky drift in a only (i.e., assuming that 
no captures occur). The result is shown in Fig. [151 If L < 1 in all the sim- 
ulations, the distribution of the asteroid semimajor axes would show a gap 
or a depression at the resonance location. On the other hand, asteroids with 
L > 1 would be the ones contributing to the total excess of asteroids inside 
the resonance. Looking at the left part of Fig. [T^] (L < 1), it is evident that 
the smaller the asteroids the larger the scaled lifetime, therefore, an excess of 
small asteroids inside the resonance with respect to the background should 
be expected. This is consistent with our analysis of the size distributions 



(Section 12.21 and Fig. [7]). It is worth mentioning that iTsiganis et al.l (120031 ) 
also found an excess of small asteroids in the 7:3 resonance with Jupiter due 
to a selective dynamical process generated by the Yarkovsky effect. 

The higher capture probability at larger eccentricities, and the diffusive 
process in eccentricity inside the resonance, explain the different eccentric- 
ity distributions observed in the resonant population compared to the non 
resonant background population. 

3.2. V-type asteroids as tracers 

V-type asteroids have unique surface colors, and they are considered to 
be members of a dynamical family formed as a result of an impact on Vesta's 
surface. There are several V-type asteroids far from the domain of the Vesta 
dynamical family in the space of proper elements, a fact that has challenged 
the understa nding of the dynamica l link between these asteroids and the 



Ve sta familv (|Nesvornv et al 



20081 ) 



Roig and Gil-HuttonI ( l2006l ). at their Table 2, provide a list of candidate 



V-type asteroids that were not identified as members of the Vesta dynamical 
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Asteroid 


ap[AU] 






H 




Aa[°] 


(oodlj Kvistaberg 


2.41^ 


^75 


0.114 


3.58 


13.2 


160 


360 


(5599) 1991 SGI 


2.41^ 


^69 


0.155 


7.23 


12.5 


175 


360 


(12027) 1997 AB5 


2.41^ 


m 


0.174 


6.22 


14.6 


183 


287 


(40733) 1999 SM17 


2.41^ 


m 


0.113 


6.80 


15.0 


167 


340 


(68064) 2000 YU67 


2.41^ 


m 


0.152 


7.13 


15.2 


179 


360 


(122441) 2000 QH126 


2.41^ 


m 


0.100 


6.92 


15.4 


177 


279 


(109721) 2001 RE54 


2.41^ 


m 


0.051 


3.93 


16.0 


217 


180 


(2798) Vergilius 


2.41^ 


m 


0.044 


5.94 


13.1 


214 


200 


(32366) 2000 QA142 


2.41^ 


m 


0.213 


6.31 


13.9 


158 


292 



Table 5: Resonant V-type ast eroids not identified a s mem bers of Vesta dynamical family. 
The first seven ones are from Roig and Gil-HuttonI ( 2006 ). Op, Cp and ip are the proper 
semimajor axes, eccentricities and inclinations, respectibely, H is the absolute magnitude, 
and (Tc and Act are the libration center and libration amplitude, respectively. 



family. From that list, we identified seven asteroids evolving in the 1:2M 
resonance, whose dynamical parameters are shown in Table [5l These aster- 
oids are very close to the region of the Vesta dynamical family, so they could 
have evolved in Cp due to the diffusive process inside the resonance, masking 
their dynamical origin as Vesta family members. The only exceptions are 
probably (3331) Kvistaberg and (109721) 2001 RE54, because their ip could 
not be explained by a diffusive process inside the resonance according to our 
numerical studies (see for example Fig. [9]). 

Looking at the photometric data from the SDSS, we identified other two 
candidate V-type asteroids among the resonant population: (32366) 2000 
QA142 and (2798) Vergilius. The orbital evolution of these two asteroids, 
together with (12027) 1997 AB5, are shown in Fig. [161 The low eccentricity 
of (2798) and the high eccentricity of (12027) are both compatible with the 
diffusive processes inside the resonance. 

The case of (32366) is a very special one, because its Cp is larger by 
approximately 0.1 with respect to the value of Cp expected for the V-type 
asteroids. Assuming a geometric albedo of ~ 0.3 — 0.4, typical of V-types, 
the diameter of asteroid (32366) is ~ 3.5 — 4.0 k m. This size is con sistent 



with a drift rate |d| 1.25 - 1.43 x 10"^ AU/Myr flRoig et al.l . 120071 ). so we 



can apply Eq. ([2]) or look at Fig. [T3]to conclude that the necessary timescale 
of resonant evolution to increase the eccentricity by 0.1 is approximately 800 
Myr. In fact. Fig. [16] shows that this asteroid has important excursions in 
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(em,^m); thus, it could have obtained its present eccentricity after several 
hundred millon years of resonant evolution. 

4. Discussion and Conclusions 

Our analysis of the size distribution of the real asteroids in the region 
of the 1:2M resonance has confirmed a significant lack of low H asteroids 
in the resonant population compared to the non resonant one. This result 
is consistent with the ratio of small to big asteroids inside and outside the 
resonance that we obtained from our simulations of fictitious asteroids, using 
different values of the drift in a induced by the Yarkovsky effect according 
to different asteroid sizes. It is worth mentioning that pure gravitational 
processes, like chaotic diffusion, do not distinguish between asteroid sizes; 
therefore in the absence of the Yarkovsky effect the size distribution in the 
resonant population should be the same as in the non resonant population. 

The diffusion in (cm, im) generates misidentification of resonant asteroids 
belonging to the Vesta and Nysa-Polana families, whose proper elements dif- 
fer from the nominal values of the families. The Massalia dynamical family 
is less affected by the resonance, probably due to the low proper inclina- 
tion of its mem bers that protects t h em fr om large diffusion in eccentricity 



and inclination. IVokrouhlicky et al.l (120061 ) modeled the orbital evolution of 



Massalia family members and found that 16% of the members that cross the 
1:2M resonance are misidentified due to the effect of the resonance. This is 
in good agreement with our result from Table |21 indicating that inside the 
resonance about 11% of the Massalia family members would not be detected 
as such. 

When the resonant asteroids reach the chaotic region at > 0.22, they 
escape from the resonance mostly by diffusion in a. A small fraction remains 
captured at higher eccentricities, but because of their large amplitude of 
cr they undergo encounters with Mars and escape from the resonance too. 
Escapes by reaching the secular resonance were also observed. For the 
small asteroids, the escape by diffusion in a is the rule, but for the big ones 
it first operates a diffusion in eccentricity, until they reach > 0.22, and 
then operates the diffusion in a. A fraction of approximately 6% of the long 
lived asteroids end by colliding with the Sun, mainly after reaching the 3:1 
resonance with Jupiter. Asteroids that tend to follow this path are the ones 
with the higher inclinations. 
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Campins et all (120101 ) discussed the origin of the NEA (101955) 1999 
RQ36, with a diameter of 580 m, and concluded that it comes from the Polana 
branch of the Nysa-Polana complex, by reaching the resonance due to the 
Yarkovsky effect over long timescales. According to our results, we suggest 
another possible path to transfer asteroids from the Polana region to a NEA 
orbit: over a small timescale, an asteroid could reach the 1:2M resonance, 
excite its eccentricity by diffusion while evolving in resonant motion, and 
finally reach the uq or the 3:1 reson a nces on a much shorter timescale than 



the one proposed by lCampins et al.l (120101 ). because for higher eccentricities 



the diffusion in semimajor axis is faster, and the secular resonance is closer 
to the Polana branch. 

The mean lifetime in the resonance of the largest members of the res- 
onant population, i.e., those not affected by the Yarkovsky effect, is quite 
long. These objects can stay evolving in the resonance over timescales of 
1 Gyr, even considering the effects of collisions with other asteroids and/or 
ejections due to close encounters with massive asteroids. This must leave 
some signature in the orbits of asteroids with diameters larger than 10 km, 
because after 1 Gyr of resonant evolution, the change in eccentricity should 
be of the order of ~ 0.04 — 0.15 (Fig. [9]). On the other hand, small resonant 
asteroids should have been captured in more recent times. 

The negligible effect of the Yarkovsky drift on a large asteroid like (142) 
Polana raises the problem of the origin of its present orbit. Ther e is an ana- 



logue case in the resonance 7:3 with Jupiter: asteroid (677) Aaltje ( iTsiganis et al. 



20031 ). One possible mechanism to explain such cases would be chaotic dif- 
fusion enhanced by close encounters with massive asteroids which, over Gyr 
timescales, might push the big asteroids into the resonance. A more efficient 
alternative, although not yet proved to be compatible with the dynamical his- 
tory of Mars, would be a slow migration of Mars orbit. This later mechanism 
of capture is dynamically possible, as we have checked with some numerical 
experiments, but it does not distinguish between asteroid sizes. 

The existence of different H distributions in the resonant and non res- 
onant populations, with an excess of small asteroids inside the resonance 
compared to the background, is an indication that the main mechanism that 
replenishes the resonance at present is the Yarkovsky effect, rather than 
either chaotic diffusion, close encounters with massive asteroids, or Mars 
orbital migration. 
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Appendix: the excess of small asteroids in the resonance 

We have defined L = trcs/T, where T is the time that would be necessary 
to cross the resonance region due to the Yarkovsky drift in a only. Taking 
a mean width for the resonance of 0.004 AU we have T = 0.004 AU/a. In 
steady state, the number of resonant asteroids, n, is 

n = N < trcs > (3) 

where 

is the injection's rate of asteroids into the resonance and < tres > the mean 
lifetime in resonance, both depending on the asteroid size. By definition of 
L we have < trcs >=< L > T for each of the four simulations. Now, taking 
into account Eqs. |3] and H] and the definition of T, we can calculate the ratio 
Us/ni between small and large asteroids inside the resonance 

_ dNs < Ls > , . 

ni ~ dNi < Li> ^ ^ 

where dNg/dNi is the ratio outside the resonance. The behavior of L as seen 
in Fig. [15] is not very different for the four simulations at the right part of 
the plot, corresponding to long lived resonant asteroids. But, looking at the 
left part of the plot (L < 1) it is evident that Lg/Li > 1 then, an excess of 
small asteroids inside the resonance with respect to outside is expected. 
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(142) Polana 



360 




20000 40000 60000 80000 100000 

time (yr) 

Figure 1: Asymmetric librations of the critical angle a = 2\~ Am — vj oi (142) Polana, 
switching its libration center at < 19.000, 30.000 and 77.000 yr. The eccentricity varia- 
tions are due to the secular forced modes. 
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Figure 2: Frequency ?;s. proper semimajor axis at the end of a numerical simulation of 
4022 real asteroids over 30,000 years. Each asteroid is represented by a dot corresponding 
to its Op and to the most representative frequency in the time evolution of a. The central 
concentration at ap ~ 2.4184 corresponds to librations and horseshoes, the straight lines 
(flp ^ 2.4174 and Up > 2.4194) correspond to circulations, and the intermediate zone 
corresponds to transitions between horseshoes and circulations. 
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Figure 3: Proper elements obtained from a numerical integration of 4022 real asteroids over 
1 Myr. Each asteroid is represented by a dot corresponding to its Up and e^. The central 
region corresponds to librations and horseshoes. The increasing width of the resonance as 
a function of e„ is evident. 
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Figure 4: Mean orbital elements, obtained from a 1 Myr running window, of 100 resonant 
asteroids from the ASTORB database numerically integrated over 1 Gyr. The concen- 
tration at a„i ~ 2.4184 AU is due to librations and horseshoes. The concentration at 
the borders of the resonance is due to the phenomenon of stickiness. The behavior at 
Cm > 0.2 is due to chaotic diffusion dirven by the overlap of several high order resonances. 
The limits of the resonance are approximately given by the shaded lines. (Low resolution 
version) 
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Figure 5: Proper elements of the population of 4022 asteroids inside and close to the 
resonance. It is possible to identify three clusters that can be associated to the families of 
Vesta, Nysa-Polana and Massalia. See Table[2]for details. The empty region over ip ^ 15° 
is due to the secular resonance fg. 
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Figure 6: Cumulative distribution of the eccentricities for the resonant and non resonant 
asteroids. The resonant asteroids have comparatively more eccentric orbits than the non 
resonant ones. 
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Figure 7: Cumulative distribution of the absolute magnitudes for two unbiased samples 
of resonant and non resonant asteroids. The two samples of 855 asteroids each, have the 
same eccentricity distribution in order to avoid possible discovery biases due to differ- 
ent perihelion distances. The resonant asteroids are comparatively smaller than the non 
resonant ones. 
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Figure 8: Time decay of the sample of 100 asteroids shown in Fig. 01 initiaUy in resonant 
motion. A power law fit gives a half-life inside the resonance of 2.5 Gyr. 
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Figure 9: Diffusion in (e™, im) of the sample of 100 asteroids shown in Fig. Diffusion is 
approximately proportional to the time span of resonant evolution. Squares indicate the 
initial values of the orbits. (Low resolution version) 
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Figure 10: Simulation including a Yarkovsky effect corresponding to a fictitious asteroid 
with a radius of some kilometers (simulation #4). The asteroid suffers a drift in semimajor 
axis up to get slicked to the resonance, then starting a diffusive process in e. The exact 
location of the resonance is indicated with a line. 
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Figure 11: The evolution of the whole population of 200 particles in the simulation with a 
Yarkovsky effect corresponding to an asteroid with a radius of some kilometers (simulation 
#4), initially located at the right side of the resonance and with a < 0. Dotted lines 
correspond to the limits of the resonance as computed in Fig. |4l (Low resolution version) 
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Figure 12: The evolution of the whole population of 200 particles in the simulation with a 
Yarkovsky effect corresponding to an asteroid with a radius of some kilometers (simulation 
^^), initially located at the left side of the resonance and with d > 0. Note the lower 
diffusion compared to Fig. [9l Particles that remain captured into the resonance over a 
short period of time, or not captured at all. do not show diffusion in (ermim)- Note also 
the diffusion in inclination for Cm > 0.22. (Low resolution version) 
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Figure 13: Total variations in Cm as a function of the total lifetime in the resonance, ^rosj 
for the simulations ^2, ^3 and #4, starting from both the left and right sides of the 
resonance. The weaker the Yarkovsky effect, the longer the time evolving in the resonance 
and the larger the total diffusion in eccentricity. The total diffusion in eccentricity follows 
an approximate power law with t^es- 
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Figure 14: Fraction of the population captured into the resonance, differentiating between 
the left and right sides populations, in the case of simulation ^^3 (i.e., with a Yarkovsky 
effect corresponding to asteroids with radii of some hundreds meters). All the four sim- 
ulations (see Table |4]) show a similar behavior: the evolution from the left and right 
populations is indistinguishable, and the tails of long lived asteroids in both cases are 
approximately the same. 
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Figure 15: Cumulative distribution of the scaled lifetime in the resonance L for all the 
simulations with Yarkovsky effect. The left and right populations appear merged. The 
lifetime in the resonance is scaled according to the time T necessary to drift by 0.004 
AU, which is an adopted mean width for the resonance. The decay of the populations is 
proportional to the Yarkovsky effect, but in the scaled time the decay is approximately 
similar for all the simulations, especially at the right part of the plot (i.e. long lived aster- 
oids). L > 1 corresponds to asteroids that have their migration slowed, on average, by the 
resonance, whereas L < 1 corresponds to asteroids that have their migration accelerated 
by the resonance. 
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Figure 16: Diffusion over 1 Gyr of three candidate V-type asteroids that evolve inside the 
1:2M resonance and were not identified as members of the Vesta dynamical family. These 
asteroids are included in the sample of 100 asteroids whose evolution is shown in Fig. El 
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